Search for: All records

An adequate characterization of internal modes of climate variability (MoV) is a prerequisite for both accurate seasonal predictions and climate change detection and attribution. Assessing the fidelity of climate models in simulating MoV is therefore essential; however, doing so is complicated by the large intrinsic variations in MoV and the limited span of the observational record. Large ensembles (LEs) provide a unique opportunity to assess model fidelity in simulating MoV and quantify intermodel contrasts. In this work, these goals are pursued in four recently produced LEs: the Energy Exascale Earth System Model (E3SM) versions 1 and 2 LEs, and the Community Earth System Model (CESM) versions 1 and 2 LEs. In general, the representation of global coupled modes is found to improve across successive E3SM and CESM versions in conjunction with the fidelity of the base state climate while the patterns of extratropical modes are well simulated across the ensembles. Various persistent shortcomings for all MoV are however identified and discussed. The results both demonstrate the successes of these recent model versions and suggest the potential for continued improvement in the representation of MoV with advances in model physics.

Modes of variability play a critical role in prediction of seasonal to decadal climate variability and detection of forced climate change, but historically many modes have been poorly simulated by coupled climate models. Using recently produced large ensembles, this work demonstrates the improved simulation of a broad range of internal modes in successive versions of the E3SM and CESM and discusses opportunities for further advances.

 

Abstract The emergence of a spatial pattern in the externally forced response (FR) of dynamic sea level (DSL) during the altimeter era has recently been demonstrated using climate models but our understanding of its initial emergence, drivers, and implications for the future is poor. Here the anthropogenic forcings of the DSL pattern are explored using the Community Earth System Model Large Ensemble (CESM-LE) and Single-Forcing Large Ensemble, a newly available set of simulations where values of individual forcing agents remain fixed at 1920 levels, allowing for an estimation of their effects. Statistically significant contributions to the DSL FR are identified for greenhouse gases (GHGs) and industrial aerosols (AERs), with particularly strong contributions resulting from AERs in the mid-twentieth century and GHGs in the late twentieth and twenty-first century. Secondary, but important, contributions are identified for biomass burning aerosols in the equatorial Atlantic Ocean in the mid-twentieth century, and for stratospheric ozone in the Southern Ocean during the late twentieth century. Key to understanding regional DSL patterns are ocean heat content and salinity anomalies, which are driven by surface heat and freshwater fluxes, ocean dynamics, and the spatial structure of seawater thermal expansivity. Potential implications for the interpretation of DSL during the satellite era and the longer records from tide gauges are suggested as a topic for future research. 

Biomass burning aerosol (BBA) emissions in the Coupled Model Intercomparison Project phase 6 (CMIP6) historical forcing fields have enhanced temporal variability during the years 1997–2014 compared to earlier periods. Recent studies document that the corresponding inhomogeneous shortwave forcing over this period can cause changes in clouds, permafrost, and soil moisture, which contribute to a net terrestrial Northern Hemisphere warming relative to earlier periods. Here, we investigate the ocean response to the hemispherically asymmetric warming, using a 100-member ensemble of the Community Earth System Model version 2 Large Ensemble forced by two different BBA emissions (CMIP6 default and temporally smoothed over 1990–2020). Differences between the two subensemble means show that ocean temperature anomalies occur during periods of high BBA variability and subsequently persist over multiple decades. In the North Atlantic, surface warming is efficiently compensated for by decreased northward oceanic heat transport due to a slowdown of the Atlantic meridional overturning circulation. In the North Pacific, surface warming is compensated for by an anomalous cross-equatorial cell (CEC) that reduces northward oceanic heat transport. The heat that converges in the South Pacific through the anomalous CEC is shunted into the subsurface and contributes to formation of long-lasting ocean temperature anomalies. The anomalous CEC is maintained through latitude-dependent contributions from narrow western boundary currents and basinwide near-surface Ekman transport. These results indicate that interannual variability in forcing fields may significantly change the background climate state over long time scales, presenting a potential uncertainty in CMIP6-class climate projections forced without interannual variability.

 

Abstract Climate models consistently project a significant drying in the Caribbean during climate change, and between 2013 and 2016 the region experienced the worst multiyear drought in the historical period. Although dynamical mechanisms have been proposed to explain drought in the Caribbean, the contributions from mass convergence and advection to precipitation minus evaporation ( P − E ) anomalies during drought are unknown. Here we analyze the dynamics of contemporaneous droughts in the Caribbean by decomposing the contributions of mass convergence and advection to P − E using observational and simulated data. We find that droughts arise from an anomalous subsidence over the southeastern Caribbean and northeastern South America. Although the contributions from mass convergence and advection vary across the region, it is mass convergence that is the main driver of drought in our study area. A similar dynamical pattern is observed in simulated droughts using the Community Earth System Model (CESM) Large Ensemble (LENS). 

The two dominant drivers of the global mean sea level (GMSL) variability at interannual timescales are steric changes due to changes in ocean heat content and barystatic changes due to the exchange of water mass between land and ocean. With Gravity Recovery and Climate Experiment (GRACE) satellites and Argo profiling floats, it has been possible to measure the relative steric and barystatic contributions to GMSL since 2004. While efforts to “close the GMSL budget” with satellite altimetry and other observing systems have been largely successful with regards to trends, the short time period covered by these records prohibits a full understanding of the drivers of interannual to decadal variability in GMSL. One particular area of focus is the link between variations in the El Niño−Southern Oscillation (ENSO) and GMSL. Recent literature disagrees on the relative importance of steric and barystatic contributions to interannual to decadal variability in GMSL. Here, we use a multivariate data analysis technique to estimate variability in barystatic and steric contributions to GMSL back to 1982. These independent estimates explain most of the observed interannual variability in satellite altimeter-measured GMSL. Both processes, which are highly correlated with ENSO variations, contribute about equally to observed interannual GMSL variability. A theoretical scaling analysis corroborates the observational results. The improved understanding of the origins of interannual variability in GMSL has important implications for our understanding of long-term trends in sea level, the hydrological cycle, and the planet’s radiation imbalance.

 

The satellite altimeter record has provided an unprecedented database for understanding sea-level rise and has recently reached a major milestone at 25 years in length. A challenge now exists in understanding its broader significance and its consequences for sea-level rise in the coming decades and beyond. A key question is whether the pattern of altimeter-era change is representative of longer-term trends driven by anthropogenic forcing. In this work, two multimember climate ensembles, the Community Earth System Model (CESM) and the Earth System Model Version 2M (ESM2M), are used to estimate patterns of forced change [also known as the forced response (FR)] and their magnitudes relative to internal variability. It is found that the spatial patterns of 1993–2018 trends in the ensembles correlate significantly with the contemporaneous FRs (0.55 ± 0.10 in the CESM and 0.61 ± 0.09 in the ESM2M) and the 1950–2100 FRs (0.43 ± 0.10 in the CESM and 0.51 ± 0.11 in the ESM2M). Unforced runs for each model show such correlations to be extremely unlikely to have arisen by chance, indicating an emergence of both the altimeter-era and long-term FRs and suggesting a similar emergence in nature. Projected patterns of the FR over the coming decades resemble those simulated during the altimeter era, suggesting a continuation of the forced pattern of change in nature in the coming decades. Notably, elevated rates of rise are projected to continue in regions that are susceptible to tropical cyclones, exacerbating associated impacts in a warming climate.

 

To provide context for observed sea level rise, the forced response (FR) in dynamic sea level (DSL) during the 20th and 21st centuries is examined in the Community Earth System Model Large Ensemble (LE). After accounting for simulation drift, which in the LE is sizable, the DSL FR is found to be complex, both in space and time. Its evolving character is suggested to arise from both the diversity and transient evolution of climate forcing agents and the slow adjustment timescales of the intermediate and deep oceans. Nonetheless, various intervals of spatially coherent change simulated for the recent past and near future are identified, and their characteristics and associated driving mechanisms are identified and discussed. The mid‐ to late‐20th century DSL FR is characterized by an hemispherically asymmetric pattern of change, with depressed rates of rise in the northern oceans. There is also a dipole of change in the Southern Ocean caused by changes in near‐surface zonal winds. Through the late 20th and early 21st centuries, a different pattern of rise emerges with elevated rates in the tropics and depressed rates at high latitudes. Zonal and interbasin variations characterize both intervals and involve the pattern and depth of anomalous ocean heat content storage and spatial contrasts in the expansion coefficient tied mainly to base state temperature, with greater rates of rise in warm regions per unit warming. The relative roles of surface flux and ocean convergence anomalies are examined.

 

Monsoon responses to eruptions over the last millennium (LM) are examined in an ensemble of climate simulations as a function of eruption hemisphere. A composite analysis reveals a particularly strong sensitivity of monsoon rainfall in the year following Northern Hemisphere (NH) extratropical eruptions. Additional analysis focusing on the 18th century eruption of Mt. Laki and idealized simulations representing an analogue Southern Hemisphere eruption (SH‐Laki) reveal monsoon responses that are approximately symmetric across hemispheres, despite exhibiting asymmetries in other aspects of the climate response. We conclude that 1) latitudinally mirrored eruptions result in approximately symmetric monsoon responses, 2) disproportionate weakening (strengthening) of NH (SH) monsoons by NH eruptions over the LM resulted in part from their relatively high latitudes, and 3) uncertainty in eruption latitude fundamentally limits our ability to accurately simulate associated monsoon and tropical precipitation responses in nature.